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Motivated by the unexplored complexity of the phase diagrams for multi-orbital Hubbard models, 
a three-band Hubbard model at integer fillings (A'^ — 4) with orbital degeneracy lifted partially by 
crystal field splitting is analyzed systematically in this work. By using single site dynamical mean- 
field theory and rotationally invariant Gutzwiller approximation, we have computed the full phase 
diagram with Coulomb interaction strength U and crystal field splitting A. We find a large region 
in the phase diagram, where an orbital selective Mott phase will be stabilized by the positive crystal 
field lifting the orbital degeneracy. Further analysis indicates that the Hund's rule coupling is 
essential for the orbital selective Mott phase and the transition toward this phase is accompanied 
by a high-spin to low-spin transition. Such a model may be relevant for the recently discovered 
pnictides and ruthenates correlated materials. 

PACS numbers: 71.10.Fd, 71.28.+d, 71.30.-|-h 



I. INTRODUCTION 

The Mott-Hubbard metal-insulator transition (MIT) 
has been a subject of great interest for decades ii Most of 
the attentions before this century have been focused on 
the one band case only, because most of the qualitative 
features of MIT have already been captured by single 
band Hubbard model, as shown by numerous studies us- 
ing the single site dynamical mean- field theory (DMFT) ^ 
Whereas in realistic materials most of the Mott transi- 
tions involve more than one band and thus exhibit multi- 
orbital features."^ Multi-orbital extension of the Hubbard 
model allows more realistic description of MIT and other 
strongly correlated physics, which contains, in general, 
much richer phase diagrams, and exotic physical phenom- 
ena. For instance, the redistribution of electrons among 
different orbitals leads to new scenarios such as orbital 
ordering;^ high-spin to low-spin (HS-LS) transition^ or 
orbital selective Mott transition (OSMT)fi etc. 

Mott transitions in multi-orbital models have been 
studied within the framework of DMFT for more than ten 
yearsi^ Previous studies show* that simple increasing 
the band degeneracy only changes the critical interaction 
strength Uc, but will not change the fundamental fea- 
tures of Mott transition, where all the degenerate bands 
undergo Mott transition simultaneously under the incre- 
ment of interaction strength. Recent studies of realistic 
materials have focused interest on the interplay between 
the MIT and orbital degeneracyi^"— A very fundamen- 
tal question raised in this field is, how the multi-orbital 
systems response to the breaking down of the orbital de- 
generacy. In such a system, it is possible that the Mott 
transitions in different orbitals happen separately, which 
is the so-called OSMT and has been suggested firstly by 
Anisimov et al.^ in the study of Ca2-3pSr2;Ru04. 

After the concept of OSMT has been proposed, lots 



of research interests have been attractedFi^Ti^i The early 
DMFT studies on this problem are focused on the two- 
band Hubbard model with half filling, which is the sim- 
plest system that may have OSMT when the band- 
widths of the two bands are different. The DMFT cal- 
culations from different groups with exact diagonaliza- 
tion (ED)— ii^ and Hirsch-Fye quantum Monte Carlo 
(HFQMC) — "= as impurity solvers converge to two es- 
sential conclusions: (1) The OSMT in two-band Hub- 
bard model is mainly induced by the bandwidth differ- 
ence which breaks the degeneracy between the two bands 
and the crystal field splitting plays a minor role here. (2) 
The emergence of OSMT is very sensitive to the feature 
of the local interaction. It is very easy to occur OSMT 
when the local interaction is rotationally invariant, while 
very difficult to happen when the local interaction breaks 
the rotational invariance. 

The OSMT in the three-band Hubbard model, 
which is more relevant to the realistic situation of 
Ca2_2;Srj;Ru04,^'^^ is not a trivial generalization of the 
two-band model. In two-band systems, the OSMT can 
only happen in half filling case, while in three-band sys- 
tems it can happen when the occupation numbers are 
2, 3 or 4. When the total occupation number is 3, the 
three-band system is half filling and actually the situa- 
tion is very similar with the two-band case. The most 
interesting case is when the occupation number is 2 or 4, 
which can be transformed between each other by particle- 
hole symmetry. Recently de' Medici et al^ proposed a 
new mechanism for OSMT, which happens in three-band 
systems with filling factor being 2 or 4. In this new sce- 
nario of OSMT, the driving force is not the difference 
of bandwidth but the crystal field splitting lifting the 
band degeneracy. And then Kita et alP'^ investigated 
how the orbital level splitting and Ising-type Hund's rule 
coupling affect the Mott transition in the case of two 



electrons per site. Their results reveal that the critical 
interaction strength separating a metallic phase and two 
kinds of insulating phases shows a non-monotonic behav- 
ior as a function of the level splitting. They suggest that 
this behavior is characteristic for 1/3 filling, in compari- 
son with the preceding results for different fillings and for 
two-band models. We note that the three-band system is 
very popular in transition metal compounds.^''' Provided 
the Fermi energy falls into the t2g bands, the tetrago- 
nal distortion will further split the t2g bands into non- 
degenerate aig band and two-fold degenerate e'g bands. 
When the occupation number is 4, the appropriate crys- 
tal field will redistribute the four electrons into aig band 
and e'g bands as 1 and 3 respectively (i.e (3,1) configu- 
ration). Thus if we neglect the correlation between aig 
and e' sub-systems, the aig band becomes a one-band 
system at half filling and e' bands become a two-band 
system at quarter filling, which will lead to an orbital se- 
lective Mott phase (OSMP) at equal bandwidth because 
it is much easier to get Mott insulator in aig band. The 
mean-field phase diagram of this model determined by 
slave-spin method^i^ contains quite a large region for 
OSMP, in which the aig band has already become Mott 
insulator while the e' bands are still metallic. 

The OSMT driven by the band degeneracy lifting is 
quite a robust phenomena determined by the interplay 
between crystal symmetry and correlation effects. In ref- 
erence [231 , the phase diagram of the three-band Hubbard 
model is mainly calculated by the slave-spin method, 
which is qualitatively correct but not accurate enough. In 
the present paper, we study systematically this t2g Hub- 
bard model with crystal field splitting by two more accu- 
rate methods: the DMFT method combined with state- 
of-the-art hybridization expansion continuous time quan- 
tum Monte Carlo (CTQMC) impurity solver^^^"^ and 
newly developed rotationally invariant Gutzwiller ap- 
proximation (RIGA) method. ^''"^■^ The metal-insulator 
phase diagram, band specific quasiparticle weight Za and 
occupation number tIq are computed by the both meth- 
ods with respect to crystal field splitting A and Coulomb 
interaction strength U. Based on these results, we mainly 
discuss three important aspects of OSMT in this system: 
(1) The crucial role of Hund's rule coupling; (2) The re- 
distribution of four electrons among aig and e' bands; (3) 
The relationship between OSMT and HS-LS spin state 
crossover. 

The rest of this paper is organized as follows: In SecHIl 
the three-band Hubbard model treated in this work is 
specified. In Sec lIII Al the main results of this paper, 
[/ — A phase diagrams for rotationally invariant interac- 
tion and SU(iV) density-density interaction are presented 
and compared with each other. In Sec lIII Bl the redistri- 
bution of electrons among different orbitals and its rela- 
tionship with the OSMT are discussed in detail. And the 
accompanying HS-LS spin state crossovers are discussed 
in Sec lIII Cl Section ITVl serves as a conclusion. 



II. MODEL 



We consider the three-band Hubbard model defined by 

H=- Y^ Uj4aaCjaa + Y. ^Ic, (1) 



E 
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where Cj^^. is an annihilation (creation) operator of an 
electron with spin a (=t, i) and orbital a (=1, 2, 3) at 
the ith site, tij is the hopping integral between site i and 
site j. The local part of Hamiltonian Hl^^ can be defined 
as follows (for the sake of simplicity, the site index i has 
been ignored): 

Hloc = - Yi^J■ - ^a)na,<, + ^ Una,\na.X 
a.a a 

+ ^ [U'Ua.anb-a + {W ~ J)na,c,nb,cy\ 

a>b,a 



a<b 



Here Ua.a = 4c 



(2) 

is the number operator, fj. is the 
chemical potential, and A^ is the energy level for orbital 
a. In the interaction terms, U{U') is the intra-orbital 
(inter-orbital) Coulomb interaction, and J is the Hund's 
rule coupling. The constrained condition U = U' + 2J 
is imposed as usual, which is valid for atomic like local 
orbitals. The above interaction terms include both the 
spin-flip and pair-hopping terms, and thus are rotational 
invariant in the spin space. In order to elucidate the Mott 
MIT in the case of four electrons per site, the chemical 
potential /i is adjusted dynamically in the simulations to 
fix the electron filling per site as A^ = 4. In the present 
paper, we focus on the i2g-like bands under the tetrago- 
nal crystal field, which are split into non-degenerate aig 
band and double degenerate e'g bands. Therefore the on- 
site energy level is assumed to be Ai y^ A2 = A3 and 
their difference is defined as A = Ai — A2. 

This lattice model (see Eq.(IT])) can be solved in the 
framework of single site DMFT— il, which neglects the 
momentum dependence of the self-energy and reduces 
the original lattice problem to the self-consistent solu- 
tion of an effective impurity model. In this paper, a 
semicircular density of states with half bandwidth D ~ 1 
is used, which corresponds to the infinite coordination 
Bethe lattice. All the orbitals have equal bandwidth and 
the energy unit is set to be 13. To solve the effective im- 
purity model, the hybridization expansion CTQMC im- 
purity solver—"— is adopted. This method allows us to 
access the strong interaction regime down to very low 
temperature. In our calculations, the system temper- 
ature is set to be T = 0.01 (corresponding to inverse 
temperature /3 = 100) unless otherwise stated. In each 
DMFT iterations, typically 4 x 10^ QMC samplings have 
been performed to reach sufficient numerical accuracy. 

In the present works, we also use the RIGA 
method^"— to crosscheck the results obtained by 



DMFT+CTQMC method. In contrast to the finite tem- 
perature DMFT+CTQMC calculations, only the zero 
temperature physical quantities can be obtained by using 
the RIGA method. Although the RIGA method can not 
access the dynamical properties of correlated systems, it 
provides a very fast and economic way to calculate the 
ground state properties, for examples, the total energy 
and occupation numbers. And it can be used as a good 
complementary method to the DMFT approach. As will 
be discussed in detail in the following sections, the phase 
diagrams and other physical properties obtained by both 
methods are in very good agreements. The implementa- 
tion details about our RIGA method have been presented 
in reference 133. 



III. RESULTS AND DISCUSSION 
A. [/ — A phase diagram 

To map out the metal-insulator phase diagram we have 
computed systematically the dependence of charge den- 
sity, quasiparticle weight, and local magnetic moment as 
the function of crystal field splitting for various Coulomb 
interaction strengths. In this subsection, we will focus 
firstly on the L'^ — A phase diagrams with both non-zero 
and zeroed Hund's rule coupling J. We will show in the 
following that the Hund's rule coupling J is extremely 
crucial for the appearance of OSMP in the three-band 
systems with four electrons. 

When the spin-flip and pair-hopping interaction terms 
are taken into full considerations, the Hamiltonian for 
three-band Hubbard model is rotationally invariant in 
the spin space, and the phase diagram contains four dif- 
ferent phases, which are metal, band insulator, OSMP 
and Mott insulator phases, respectively. The obtained 
L/ — A phase diagrams for J — U/A case are illustrated 
in Figdl The upper panel is the results obtained by 
DMFT-fCTQMC method, and the lower panel is ob- 
tained by RIGA method. From the above figure, we can 
easily see that the consistency between both two methods 
is quite excellent, except the phase boundary obtained 
by RIGA method is a bit higher, which is similar with 
the situations found in one and two-band models . ^^ ' ?? 
Comparing with the similar phase diagrams obtained by 
slave-spin method in the literatures ^^i^ there are two 
important differences. Firstly, the critical Uc for Mott 
transition at A = is much lower in the phase diagram 
obtained by DMFT+CTQMC method. Secondly, in the 
phase diagram obtained by both DMFT-fCTQMC and 
RIGA methods, the phase boundary between the metal 
phase and band insulator phase depends on the interac- 
tion strength U monotonically, however in the phase di- 
agram obtained by slave-spin method, it decreases firstly 
and then increases i^ This difference is mainly due to the 
over simplified treatment of Hund's rule coupling terms 
in the slave-spin methodi^ 

The general shape of the phase diagram can be eas- 
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FIG. 1. (Color online) Calculated phase diagram of three- 
band Hubbard model with rotationally invariant interactions 
in the plane of Coulomb interaction U {J = f//4) and crystal 
field splitting A (A = Ai — A2). Upper panel: Calculated 
by DMFT-I-CTQMC method at finite temperature T = 0.01. 
Lower panel: Calculated by RIGA method at zero temper- 
ature. In all the calculations, the chemical potential /i is 
adjusted dynamically to fulfill the total occupation number 
condition (Af = 4). The label "(2,2) Metal + Band Insulator" 
means the two-fold degenerate bands (e^ states) are metallic 
while the non-degenerate band {aig state) is insulating, and 
(2,2) means corresponding orbital occupancies. All the other 
labels have similar explanations. The pink zone means LS 
state and the orange zone means HS state respectively. 



ily understood by considering two limiting cases: (1) For 
A = 0, the model reduces to a fully degenerate three- 
band Hubbard model with total filling iV = 4, which 
undergoes a Mott transition around U = 10.0; (2) For 
the uncorrelated limit {U ~ 0), when A > 0, a simple 
transition from metal to band insulator can be observed 
at A = 2.0 with the fully occupied e' bands and empty 
aig band. Another transition can also be seen on the neg- 
ative side at A = —1.0, after which the non-degenerate 
aig band is fully occupied band insulator and the double 
degenerate e' bands are metallic at half filling. 

When both the Coulomb interaction U and crystal field 
A are non-zero, the phase diagram is quite complicated 
and very different for A > and A < 0. Now let's make 
a further discussion about this phase diagram. If the 
crystal field splitting A > 0, the isolate aig band is lifted 
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FIG. 2. (Color online) Quasiparticle weights Za as a function 
of Coulomb interaction U for selected crystal field splitting 
A — 1.0. The calculations are done by the RIGA method at 
zero temperature. 



up and the double degenerate e' bands are pushed down. 
This part of phase diagram can be divided vertically into 
two regions: (1) When 2.0 > A > 0.0, the system under- 
goes successive transitions from metal to Mott insulator 
through an OSMP as the Coulomb interaction strength 
increases; (2) When A > 2.0, the system undergoes a dif- 
ferent type of successive transitions from band insulator 
to Mott insulator through metal and OSMP phases in 
sequence. On the other hand, if the crystal field splitting 
A < 0, i.e, the oig band is lower than the two-fold degen- 
erate e'g bands, this part of phase diagram can be divided 
into three different regions: (1) For —1.0 < A < 0.0, with 
weak Coulomb interaction both bands have fractional fill- 
ing and are metallic. As the increment of Coulomb inter- 
action U, the system undergoes a transition from a fully 
metallic phase to an insulator phase, after which the oc- 
cupation numbers for both aig and e' bands are two. 
Since the aig band has no degeneracy, it is completely 
occupied and becomes band insulator in this situation. 
While at the same time the e' bands with double degen- 
eracy become Mott insulator with half filling. It is called 
(2, 2) Mott insulator phase throughout this paper. (2) 
When —1.4 < A < —1.0, with weak Coulomb interaction 
the aig band is already fully occupied and the double 
degenerate e'g bands are half-filled and metallic. With 
the increment of Coulomb interaction {/, the system first 
becomes fully metallic phase with fractional filling fac- 
tors for all the bands, and then goes back to the original 
phase after a "reentrance" transition. Finally it becomes 
the (2, 2) Mott insulator phase for U > 2.0. (3) When 
A < —1.4, the aig band remains fully occupied regard- 
less of the interaction strength and the system reduces 
to a two-band system with half filling, which undergoes 
a typical Mott MIT around U = 2.0 as determined by 
DMFT-f-CTQMC method. 

We note that the effect of Coulomb interaction U have 
two important effects for multi-orbital systems. One is 
to reduce the quasiparticle weight, the other one is to 
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FIG. 3. (Color online) Calculated phase diagram of three- 
band Hubbard model with SU(A'^) density-density interaction 
in the plane of Coulomb interaction U {J = 0.0) and crystal 
field splitting A (A = Ai — A2). Upper panel: Calculated 
by DMFT-I-CTQMC method at finite temperature T = 0.01; 
Lower panel: Calculated by RIGA method at zero tempera- 
ture. In the calculations, the chemical potential /i is adjusted 
dynamically to fulfill the total occupation number condition 
(TV = 4). The label "(2,2) Metal -I- Band Insulator" means 
the two-fold degenerate bands {e'g states) are metallic while 
the non-degenerate band (aig state) is insulating, and (2,2) 
means corresponding orbital occupancies. All the other labels 
have similar explanations. The pink zone means LS state and 
the orange zone means HS state respectively. 



redistribute the electrons among different bands. The 
interplay between these two effects determines the main 
structure of the above phase diagram. Next we will focus 
on the first effect, and the second effect is less emphasized 
before and will be discussed in detail in the following 
subsections. 

The band specific quasiparticle weights Za as a func- 
tion of Coulomb interaction U are shown in Figj^l The 
crystal field splitting is fixed to A = 1.0. For simplic- 
ity, only the results obtained by RIGA method, which 
are consistent with those obtained by DMFT-f CTQMC 
method, are displayed in this figure. It is apparent that 
the quasiparticle weights Za decrease monotonously from 
1.0 to 0.0 when Coulomb interaction strength increases. 
As [/ < 3.0, the Za are large than 0.1, and both the 
e' and aig bands are metallic. As U > 3.0, Zi turns 



to zero while Z2{— Z3) is still finite. It means that the 
aig band undergoes a Mott MIT and turns into insula- 
tor at U = 3.0, and at the same time the e' bands keep 
metallic, which is the so-called (3,1) OSMP in the phase 
diagram. When Coulomb interaction strength continues 
to increase {U > 11.0 for the RIGA method and U > 9.0 
for the DMFT-I-CTQMC method), the e^ bands undergo 
another Mott MIT and the system becomes the (3,1) 
Mott insulator phase finally, then all Za will approach 
zero. 

The Hund's rule coupling J has enormous infiuence on 
the metal-insulator phase diagram for multi-orbital Hub- 
bard modeli^i^ It is noted that if we neglect the spin- flip 
and pair-hopping terms, only keep the Ising-type Hund's 
rule coupling term J^ (= J), we can obtain very simi- 
lar phase diagram, albeit the phase boundary is shifted 
downwards slightly. It is apparent that the Ising-type 
Hund's rule coupling term J^ is the key requirement for 
OSMT and such a rich phase diagram. In order to reveal 
the underlying physics more clearly we also did calcula- 
tions for the phase diagram of similar three-band Hub- 
bard model without Hund's rule coupling terms (J — 0). 
Thus we can compare the [/ — A phase diagrams for the 
case with finite J (see FiglT]) and SU(A^) J = case (see 
Figini). Obviously the latter is much simpler and has no 
OSMP in the whole phase diagram. 

The calculated U — A phase diagram for J = case 
is illustrated in Fig|3l The upper panel is the calculated 
results obtained by DMFT+CTQMC method, and the 
lower panel is obtained by RIGA method. These two 
methods give almost identical results again. Firstly we 
concentrate on the non-interaction case {U — 0). When 
the crystal field splitting is positive (A > 0), an obvious 
transition can be seen at A = 2.0 after which the system 
becomes band insulator. On the other hand, when crys- 
tal field splitting is negative (A < 0), similar MIT can 
be found at A — —1.0. When A > 0, the phase diagram 
only consists of metal and band insulator phases. The 
lower left region is metallic, and the upper right region 
is band insulator phase. The OSMP is disappeared com- 
pletely. When A < 0, the phase diagram contains metal, 
band insulator, and Mott insulator phases, but the char- 
acteristic "tip" which can be clearly seen at U ^ 1.0 
and A ~ —1.2 in Fig [T] vanishes, and the phase boundary 
between metal plus band insulator and (2,2) Mott insula- 
tor is shifted upward significantly. Thus, it is concluded 
that the Hund's rule coupling terms J has played a key 
role on the phase diagram and finite Jz is the minimal 
requirement to drive an OSMT. 



B. Redistribution of electrons 

In this subsection, we will discuss the second important 
effect of Coulomb interaction U on multi-orbital systems 
with broken symmetry: the redistribution of electrons 
among different orbitals. In order to understand the in- 
triguing physics contained in the U — A phase diagrams 
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FIG. 4. (Color online) Orbital filling as a function of Coulomb 
interaction U for selected crystal field splitting A values (rota- 
tionally invariant interaction case) . The results are calculated 
by DMFT-I-CTQMC method at finite temperature T = 0.01. 
In this figure, only the electron densities of non-degenerate 
band (aig band) are shown. The arrows correspond to phase 
transition points. 
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FIG. 5. (Color online) Orbital filling as a function of crystal 
field splitting A for selected Coulomb interaction U values, 
(a) U = 1.0; (b) U = 3.0; (c) U = 7.0; (d) U = 3.0. In (a)-(c), 
the interaction term is rotationally invariant with non-zero 
spin-flip and pair-hopping terms, whereas in (d) the interac- 
tion term is SU(A'') scheme with J = 0. All the calculations 
are done by DMFT+CTQMC method at finite temperature 
T — 0.01. In this figure, only the electron densities of non- 
degenerate aig band are shown. The arrows correspond to 
phase transition points. 



more clearly, we further plot the occupation numbers of 
non-degenerate aig band as functions of Coulomb inter- 
action U and crystal field splitting A in FigHHHl respec- 
tively. 

In FiglU the evolution of the orbital occupancy with 
the Coulomb interaction strength is shown. The vertical 
arrows denote phase transition points. In this figure, only 
the orbital filling of the non-degenerate band {aig band) 
is plotted. For negative crystal field, aig band is much 
lower in energy, and the effect of Coulomb interaction de- 
pends on the value of the crystal field A in the following 
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FIG. 6. (Color online) Orbital filling as a function of 
Coulomb interaction U for selected crystal field splitting val- 
ues (SU(Af) J = case). All the calculations are done by 
DMFT+CTQMC method at finite temperature T = 0.01. In 
this figure, only the electron densities of non-degenerate aig 
band are shown. The arrows correspond to phase transition 
points. 



way. For > A > —1.0, the effect of correlation effect is 
to transfer electrons from e' bands to the aig band until 
it is fully occupied and becomes band insulator, as shown 
by the blue curve in FigH) For -1.0 > A > -1.4, the aig 
band is already fully occupied in the non-interacting case 
and the effect of the correlation is non-monotonic. As 
the increment of the repulsive interaction U and Hund's 
rule coupling J, the occupation number of aig band first 
drops due to the Hund's rule coupling and then returns 
back to be fully occupied. We note that the interesting 
non-monotonic behavior of the occupation is the conse- 
quence of the interplay between Hund's rule coupling J, 
which favors even distribution of the electrons, and the 
Coulomb repulsive interaction U, which intends to in- 
crease the occupation difference between orbitals. There- 
fore, this behavior disappears when the Hund's rule cou- 
pling has been set to be zero, as shown in figure [51 While 
for A < —1.4, the charge distribution will not be affected 
by the correlation effect and the occupation number of 
aig band keeps constant as the increment of U . 

While the situation is very different for A > 0, in this 
case the energy level of aig band is higher and the oc- 
cupation number is less than 0.5 in the non-interacting 
case. As we can see from FiglH the increment of the in- 
teraction strength will pump the electrons from e' bands 
to aig band again until it reaches half filling and be- 
comes Mott insulator. When the crystal field strength 
A is smaller than 2.0 all the bands are metallic, and the 
system becomes OSMP when the aig band reaches half 
filling. While for A is larger than 2.0, the system starts 
from a typical band insulator with fully occupied e' and 
empty aig bands. With the increment of Coulomb inter- 
action, it first becomes metallic when the aig band get 
partially populated and finally goes into the OSMP. 

In FigO (a)-(c), we plot the occupancy of aig band 
under fixed Coulomb interaction strength C/ as a func- 



FIG. 7. (Color online) The calculated effective local magnetic 
moment v(S|) for three-band Hubbard model with rotation- 
ally invariant interaction terms. The calculations are done by 
DMFT-I-CTQMC method at finite temperature T = 0.01. In 
this figure, only the calculated results for representative crys- 
tal field splitting values (A =-2.5, 1.0, and 2.5) are presented. 
The pink region denotes LS phase, while orange region de- 
notes HS phase. 



tion of crystal field splitting A. For weak Coulomb in- 
teraction {U = 1.0), when the crystal field splitting is 
increased from very negative to very positive value, two 
plateaus can be found in the occupancy of aig band, 
which correspond to fully occupied and empty situations 
respectively. The occupation numbers decrease smoothly 
between the two plateaus corresponding to the metallic 
phase. For intermediate interaction strength with non- 
zero Hund's rule coupling J, there is additional plateau 
with occupancy being half filling, which corresponds to 
the OSMT phase which is completely absent when J = 0. 
The redistribution of the electrons among the three 
orbitals is the key point for the OSMT in this system and 
can be understood by the effect of Hund's rule coupling, 
which favors the HS state with the (3,1) configuration. 
As for SU(A'') J = case, the redistribution of electrons 
is much more easy to be understood. As shown in FiglHl 
the correlation effect induced by the Coulomb interaction 
increases the occupation of aig band for A < case and 
decreases it for A > case, which is consistent with the 
results obtained by Hartree-Fock mean-field method. 



C. HS-LS transition 

In this subsection, we focus on the magnetic proper- 
ties of the three-band model during the phase transitions 
which have been reported rarely in the references. In 
Fig|7] we plot the evolution of the mean instantaneous 
moment, which is defined as ^J {SI), as the function 
of Coulomb interaction strength U with selected crystal 
field splitting A. We find that although the overall be- 
havior looks quite similar for positive and negative crystal 
field splitting, the underlying physics are very different 



for the two cases. 

For the negative crystal field splitting A < 0, with the 
increment of Coulomb interaction strength and Hund's 
rule coupling, the system is approaching the (2, 2) con- 
figuration (two electrons in both aig and e' bands). The 
HS state can be reached smoothly when the two elec- 
trons in the e' bands fall into the triplet state due to 
the Hund's rule coupling. At the mean while, the Mott 
transition happens when the spin state gets close to the 
HS state. 

Nevertheless, the situation on the positive crystal field 
splitting side A > is very different. Under weak 
Coulomb interaction, the occupation numbers of the sys- 
tem are more close to the (4, 0) configuration (all the 
four electrons are in e' bands). With the increment of 
Coulomb interaction strength and Hund's rule coupling, 
when the crystal field splitting is still not strong enough 
to drive the system into the band insulator, the system 
undergoes two sequential transition. The first transition 
is more related to the charge transfer from the e' bands 
to aig band and the system changes from (4, 0) to (3, 1) 
configuration afterwards. This transition comes from the 
competition between the crystal field splitting, which fa- 
vors the (4, 0) configuration, and the Hund's rule cou- 
pling, which favors the (3,1) configuration. After the first 
transition, which manifests itself in FiglT] as a clear kink 
in the \/{Sf} — U curve, the aig band is half filling and 
the e' bands are only quarter filling. Since the system 
with half filling is always much closer to a Mott transi- 
tion, the aig band becomes Mott insulator first after the 
first transition, while the e' bands still remain metallic. 
A second transition happens by further increasing the 
Coulomb interaction together with the Hund's rule cou- 
pling, after which the system becomes Mott insulator for 
all the bands with the spin state being very close to a 
pure HS state. 



level splitting and the Hund's rule coupling affect the 
Mott transition in a system with four electrons per site. 
We obtain the following conclusions. 

Firstly, the Hund's rule coupling J (more specifically, 
the Ising-type Hund's rule coupling J^) is the minimal 
requirement to induce OSMT and stabilize OSMP. In 
the phase diagram for three-band Hubbard model with 
rotationally invariant interactions, the OSMP covers a 
wide parameter range. While in the phase diagram for 
three-band Hubbard model without Hund's rule cou- 
pling terms, though the Mott MIT occurs when the 
Coulomb interaction strength reaches critical values, but 
the OSMP is absent. 

Secondly, the interplay between crystal field splitting 
and Hund's rule coupling terms induces the redistribu- 
tion of electrons among the three bands and leads to the 
complex phase diagrams. For examples, the (4,0) elec- 
tronic configuration is corresponding to a band insulator 
phase, (3,1) configuration leads to OSMP or Mott insu- 
lator phase, and (2,2) configuration leads to metal plus 
band insulator or Mott insulator plus band insulator. 

Thirdly, the appearance of the OSMP in this system is 
always accompanied by a HS-LS spin state crossover. In 
the OSMP, the electron distribution keeps in (3,1) con- 
figuration, which strongly favors the HS state and lowers 
the Hund's rule energy. When the Coulomb interaction 
is reduced, the OSMP collapses to metal phase accompa- 
nied by a HS-LS transition. 

Finally, three-band models arise in a number of other 
physically important contexts, including doped Cgo, 
ruthenates, etc. Extending our results to models with 
more realistic density of states is of high priority for fu- 
ture research. 
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